Spontaneous spin polarized tunneling current through a quantum dot array 
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We show theoretically that a strongly spin-polarized 
current can be generated in semiconductors by taking 
advantage of the ferromagnetic phase of a quantum dot 
array (QDA). A Hubbard model with coupling to leads is 
used to study the tunneling current of the QDA system 
as a function of gate voltage. Due to the weak interdot 
coupling and strong Coulomb repulsion, it is found that 
a ferromagnetic phase exists in QDA within a window of 
gate voltage. Therefore QDA can be used as a spin filter 
to detect and control spin states in quantum information 
devices. 

Recently, spintronics and quantum information 
(QI) processing have attracted a great deal of 
attention. [1] A good QI system should provide 
well denned quantum computational space, pre- 
cise quantum-state preparation, coherent quantum 
manipulation, and state detection[2]. Solid state 
devices based on modern advanced semiconductor 
techniques opened up the possibility of fabricating 
large integrated networks which would be required 
for the realization of quantum computation [3, 4]. 
Both charged states [5] and spin states [6] of elec- 
trons have been proposed to carry the quantum in- 
formation. Because the decoherence time of spin is 
much longer than that of charge [6], using the spin 
is more promising for quantum information process- 
ing. Kane has proposed to use the nuclear spin as 
the quantum bit, since its decoherence time is much 
longer than that of the electron spin[7]. Neverthe- 
less, its manipulation is difficult due to the weak 
coupling between nuclear spin and electron spin (hy- 
perfine interaction). Thus, the use of spin states of 
electrons in semiconductors remains a viable option. 
One of the challenging problems is the preparation 
of electrons in semiconductors with well defined spin 
state. 

Controlling electron spin states, such as coherent 
manipulation and filtering, becomes crucial in the 
implementation of quantum computer. DiVincenzo 
[8] has suggested to use the spin filter effect to ma- 
nipulate spin states. Semiconductor quantum dots 
with local magnetic field can be used as spin filter 
and momory (read out) device [9]. Ferromagnetic 
semiconductor materials can also be used as spin fil- 
ters[10,ll]. EuO and EuS have been suggested as 
spin filters by DiVincenzo [8], but the compatibitity 
of these magnetic materials with conventional semi- 
conductors like GaAs or Si or Ge is unclear. The 
spin polarization in tunneling current has exceeded 



99% for EuO and EuS, and up to 90% in BeMnZnSc 
[11]. 

Here, we propose to use a narrow band QD ar- 
ray (QDA) weakly coupled to leads as a spin filter. 
We find that electrons injected from the leads into 
the QDA within a small window of the applied gate 
voltage will favor a ferromagnetic state, as a result 
of the strong electron correlation. In other words, 
the QDA functions as a spin filter, which may be 
used to detect and manipulate spin states. We also 
expect that the QDA has potential applications in 
spintronics. 

The device is described by the Hamiltonian, 

H = T,k t a e ^ a i t a a K^ + Ep^Ep^A,^ + 

Et,k, CT ^.kfl^o-^^ + Ei,p,<7 V i,P b p,a d i,v + h.C. + 

Ei,cr E o n ^ + EiJ *i,i4<A ff + Ei, CT U^m-v, 
where the first two terms describe the left lead 
and right lead, respectively. The third and fourth 
terms describe the coupling between the quantum 
dot (QD) and the two leads. The fifth and sixth 
terms describe the energy level of the quantum dot 
and interdot coupling. The last term describes the 
intradot Coulomb interaction. We take into account 
only one energy level for each dot and the nearest- 
neighbor coupling between dots (Uj = —t for near- 
est -neightbor For small size QDs, the energy 
difference between the ground state and the first ex- 
cited state is much larger than t and the Coulomb 
interaction U . Therefore, it is a good approximation 
to consider just one energy level in each dot. 

Using Keldysh's Green function method[12] , we 
obtain the spin-dependent tunneling current 



J* = ~ ^[fL-fRW^ImG^,^), (1) 



where lo^ denotes the energy of the electron in the 
leads with wave vector k. /l = /(wk — £tt) an d 
fn. = /( w k — Hr) are the Fermi distribution func- 
tion of the left lead and right lead, respectively, hl 
and fin are the chemical potentials in the left and 
right leads, respectively. They are related to the ap- 
plied bias, V a by /i£ — = eV a . For simplicity, 
we assume that the QD couples with the left and 
right leads symmetrically, although it is straight- 
forward to extend to the case with asymmetric cou- 
pling. Vk = Ej ^j,k e_ ^ wn 6 r 6 k|| is the pro- 
jection of electron wave vector k in the QDA plane 
and Rj is the position of the j-th QD. We propose 
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a setup in which a small bias V a (eV a is comparable 
to t) is applied, and scan the gate voltage V g (which 
serves the purpose of tuning the QD energy level E 
relative to the Fermi level in the leads) in order to 
observe the spin-dependent current. 

The calculation of tunneling current is entirely 
determined by the retarded Green function for the 
QDA. Finding the spin-polarized Green function for 
the Hubbard model was considered a challenging 
problem [13]. It is well known that the retarded 
Green function obtained within the Hubbard ap- 
proximation does not support any magnetic ordcr- 
ing[14]. However, Harris and Lange[15] showed that 
ferromagnetic ordering can be a stable state in the 
3D Hubbard model by introducing the spin depen- 
dent band shift. This mechanism plays the crucial 
role for determining the ferromagnetic state[15-18]. 
Bcencn and Edwards[13] used the approach devel- 
oped by Roth[16] to study the 2D Hubbard model for 
the normal and superconducting state of Cu02- The 
quasi-particle excitation energy they obtained is in 
very good agreement with quantum Monte Carlo cal- 
culations [19]. Therefore we adopt Roth's procedure 
to calculate the retarded Green function Gk M ,o-(w), 
while taking into account the coupling between QDA 
and leads. In the weak-coupling limit (t <C U), we 
obtain 

k||,<Tla;j ~~ uj-E ~ e fe (1 - n_ CT ) - n- a W^,,,- a + iF(k) 



+ 



w - E - U - e k n- a - (1 - n_ CT )VK k a + iT{k) ' 



where r(k, u>) denotes the the tunneling rate from 
the QDA to the leads . It is cumbersome to fully 
include the tunneling rates as a function of the 
wave vector and bias. We treat I as a constant 
parameter, even though it can be determined via 
a numerical method [20]. This approximation is 
valid for the small range of applied bais, because 
the potential barrier between leads and QDA is 
high. In the Coulomb blockade regime, it is ade- 
quate to consider the coupling between the QDA 
and the leads within the Hatree-Fock approxima- 
tion[21] (leading order). Here we consider a square 
lattice with lattice constant a. The energy dis- 
persion of electrons in the QDA is then given by 
e(k||) = —2t[cos(k x a) + cos(k y a)]. k|| is restricted in 
the first Brillouin zone of the 2D lattice. Wkn,-<7 de- 
notes the spin-dependent band shift, which is given 
by n a (l -n a )Wu n ,„ = w 0t<T + Wi ;Cr e(k||), where w ^ 
denotes the electron hopping correlation, while u>i i<T 
consists of three terms, which represent the density 
correlation, spin correlation, and spin-flip correla- 
tion, respectively [16]. Herrmann and Nolting [22] 
have proved that the effect due to w\, a is small for 
the ferromagnetic state of a body central cubic lat- 
tice. This implies that the electron hopping cor- 
rection can maintain the stability of ferromagnetic 



state in 3D. Nevertheless, is kept in the present 
calculation in order to obtain more accurate result. 

When the chemical potential in the left lead is 
lower than E + U and U 3> t, the effect due to 
the high energy pole of the Green function can be 
ignored. Consequently, the infinite U limit can be 
used in the calculation of Wk^-o- [16] which yields 



-4tni, CT _ , nf >tT (l - 7y) + ni, CT ni,_ CT 



with 



ni, 



-1 V e(k) f 
2tt ^ 4t / 



du>[f L (u) + f R (w)]ImGk..< r (w). 



The number of electron per dot is calculated by 
n ° = [ Mh(u) + /«M]/mG lr H. (2) 
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At zero temperature, the integral over w can be 
carried out analytically. Thus, we can obtain n a and 
Th — (j by solving two coupled one-dimensional integral 
equations self-consistently. The parameters used in 
our calculations are E = eV g + {^l + Mij)/2 + At 
(without applied bias), U — 20i, a — 200 A, and the 
effective mass of electrons in leads (assumed to be 
GaAs) to* = 0.067m e . 

The spin-dependent electron occupancy n a {a =| 
,|) as a function of the total electron occupancy 
n = nj + nj, at zero temperature for three different 
applied voltages is shown in Fig. 1 ; solid line denotes 
eV a = O.lt, dashed line denotes eV a — OM, and dot- 
ted line denotes eV a — t. We see a bifurcation for 
n T and n L at n = 0.364,0.367,0.373 for eV a = O.lt, 
0M, and t, respectively, where the system becomes 
spin polarized. The system remains ferromagnetic 
for n within a small window and reverts to the para- 
magnetic state at n w 0.4, beyond which the pseudo- 
equilibrium condition[16] can no longer be satisfied 
with n-f 7^ n j. As the applied bias is increased, the 
spin polarization and the domain that maintains the 
spin polarization are reduced. Eventually, the ap- 
plied bias will totally destroy the spin polarization 
of the system. 

Fig. 2 shows the spin-dependent tunneling cur- 
rent as a function of gate voltage (V g ) for various 
strengths of applied bias. We define the spin polar- 
ization of the current as V s = (Ji — Jj,)/(J-f + Jj,). 
We see that the maximun V s = 0.4816, 0.3732, and 
0.3197 for eV a = 0.1,0.5* and t. This is much bet- 
ter than the value (less than 1%) achieved by us- 
ing ferromagnetic metals in contact with semicon- 
ductors, because the conductivity of metal is much 
higher than that of semiconductor. [23,24] Although 
using magnetic semiconductor to replace the ferro- 
magnetic metal, the best value of V s achieved is near 
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90% [14], it remains to be seen if similar idea can be 
applied to III-V and group- VI semiconductors. The 
spin dependent tunneling current shown in Fig. 2 
implies that we can readily manipulate the spin po- 
larization of the tunneling current by the gate volt- 
age without introducing magnetic field or magnetic 
dopants. Although the value of the spin polariza- 
tion (P s ) obtained here is not very high, we believe 
that P s can be enhanced by using coupled multiple 
layers of QDAs, since in the 3D Hubbard model the 
ferromagnetic phase is stable over a wider range of 
n and nj (for majority carrier) can approach 1 . [16] 
According to Eq. (1), the tunneling current J± (for 
minority carrier) can be reduced to zero as a result 
of the factor (1 — n-f), and P s can approach 1. 
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Figure Captions 

Fig. 1: Spin-dependent electron occupancy versus 
total occupancy n for various strengths of applied 
bias; solid line (eV a = O.lt), dashed line (eV a = 
0.5i), and dotted line (eV a = i). 

Fig. 2: Spin dependent tunneling current as a 
function of gate volatge (V g /4t) for various strengths 
of applied bias; solid line (eV a = O.lt), dashed line 
(eV a — 0.5i), and dotted line (eV a = t). 
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